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By improving resolution of home-made mutual-inductance measurements technique, a pronounced step-like 
structure (with the number of steps n — 4 for all AC fields) has been observed in the temperature dependence 
of AC susceptibility in artificially prepared two-dimensional Josephson Junction Arrays (2D-JJA) of unshunted 
Nb — AlOx — Nb junctions with /3l(4.2A') = 30. Using a single-plaquette approximation of the overdamped 
2D-JJA model, we were able to successfully fit our data assuming that steps are related to the geometric 
properties of the plaquette. The number of steps n corresponds to the number of flux quanta that can be 
screened by the maximum critical current of the junctions. The steps are predicted to manifest themselves in 
arrays with the inductance related parameter Pl{T) matching a "quantization" condition /3l(0) = 27r(n + 1). 



PACS: 74.25.Ha, 74.50.+r, 74.80.Bj 

1. Introduction. Many unusual and still not com- 
pletely understood magnetic properties of Josephson 
Junction Arrays (JJAs) continue to attract attention of 
both theoreticians and experimentalists alike (for recent 
reviews on the subject see, e.g. [Q |3 El ED and further 
references therein). In particular, among the numerous 
spectacular phenomena recently discussed and observed 
in JJAs we would like to mention the dynamic temper- 
ature reentrance of AC susceptibility 2 (closely related 
to paramagnetic Meissner effect 0]) and avalanche-like 
magnetic field behavior of magnetization 0] [5] (closely 
related to self-organized criticality (SOC) [HIED- More 
specifically, using highly sensitive SQUID magnetome- 
ter, magnetic field jumps in the magnetization curves 
associated with the entry and exit of avalanches of tens 
and hundreds of fluxons were clearly seen in SIS-type 
arrays 5 . Besides, it was shown that the probabil- 
ity distribution of these processes is in good agreement 
with the SOC theory 0. An avalanche character of 
flux motion was observed at temperatures at which the 
size of the fluxons did not exceed the size of the cell, 
that is, for discrete vortices. On the other hand, us- 
ing a similar technique, magnetic flux avalanches were 
not observed in SNS-type proximity arrays |S] despite a 
sufficiently high value of the inductance L related crit- 
ical parameter (3l = 2ttLIc/$o needed to satisfy the 
observability conditions of SOC. Instead, the observed 
quasi-hydrodynamic flux motion in the array was ex- 
plained by the considerable viscosity characterizing the 
vortex motion through the Josephson junctions. 



In this paper we present experimental evidence for 
manifestation of novel geometric effects in magnetic 
response of high-quality ordered 2D-JJA. By increas- 
ing the resolution of our home-made mutual-inductance 
measurements technique, we were able to observe for the 
first time a fine, step-like structure (with the number of 
steps n — 4 for all AC fields used in our experiments) in 
the temperature behavior of AC susceptibility in artifi- 
cially prepared 2D-JJA of unshunted Nb — AlO x — Nb 
junctions. Using a single-plaquette approximation of 
the overdamped 2D-JJA model, we show that the num- 
ber of steps n corresponds to the number of flux quanta 
that can be screened by the maximum critical cur- 
rent of the junctions and as a result steps will man- 
ifest themselves in arrays with the inductance related 
parameter /3l(T) matching a "quantization" condition 
/? L (0) = 27r(n + l). 

2. Experimental results. To measure the com- 
plex AC susceptibility in our arrays with high preci- 
sion, we used a home-made susceptometer based on the 
so-called screening method in the reflection configura- 
tion |Ul[roi[ri]. The experimental system was calibrated 
by using a high-quality niobium thin film. Previously, 
we have shown that the calibrated output of the complex 
voltage in this experimental setup corresponds to the 
true complex AC susceptibility (for more details on the 
experimental technique and setups used in our study, 
see 

Measurements were performed as a function of the 
temperature T (for 1.5K < T < 15K), and the am- 
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plitude of the excitation field h ac (for ImOe < h ac < 
lOOe) normal to the plane of the array. The frequency 
of AC field in the experiments reported here was fixed 
at 20kHz. The used in the present study unshuntcd 
2D-JJAs are formed by loops of niobium islands (with 
T c = 9.25K) linked through Nb - AlO x - Nb Joseph- 
son junctions and consist of 100 x 150 tunnel junctions. 
The unit cell has square geometry with lattice spac- 
ing a = 46/iTO and a single junction area of 5 x 5/xm 2 . 
Since the inductance of each loop is L = i^a = QApH 
and the critical current of each junction is Ic(4k.2K) = 
150/j,A, we have (3l{^-2K) = 30. Recall that parameter 
(3l(T) = 2nLIc(T)/$ Q (where $ is the magnetic flux 
quantum) is proportional to the number of flux quanta 
that can be screened by the maximum critical current 
in the junctions. 

It is important to mention that magnetic field depen- 
dence of the critical current of the array (taken at T = 
4.2K) on DC magnetic field Hd c (parallel to the plane 
of the sample) exhibited E] a sharp Fraunhofer-like 
pattern characteristic of a single-junction response, thus 
proving a rather strong coherence within arrays (with 
negligible distribution of critical currents and sizes of 
the individual junctions) and hence the high quality of 
our samples. 

The observed temperature dependence of the real 
part of AC susceptibility for different AC fields is shown 
in Fig. H A pronounced step-like structure is clearly 
seen at higher temperatures. The number of steps n 
does not depend on AC field amplitude and is equal 
to n = 4. As expected 011111121, for h ac > AOmOe 
(when the array is in the mixed-like state with practi- 
cally homogeneous flux distribution) the steps are ac- 
companied by the previously observed reentrant behav- 
ior with x'C^j h ac ) starting to increase at low tempera- 
tures. 

3. Discussion. To understand the step-like be- 
havior of the AC susceptibility observed in unshuntcd 
2D-JJAs, in principle one would need to analyze in de- 
tail the flux dynamics in these arrays. However, as we 
have previously reported 0E1E1, because of the well- 
defined periodic structure of our arrays with no visible 
distribution of junction sizes and critical currents, it is 
quite reasonable to assume that the experimental results 
obtained from the magnetic properties of our 2D-JJAs 
could be understood by analyzing the dynamics of just 
a single unit cell (plaquette) of the array. As we shall 
see, theoretical interpretation of the presented here ex- 
perimental results based on single-loop approximation, 
is in excellent agreement with the observed behavior. 
In our analytical calculations, the unit cell is a loop 
containing four identical Josephson junctions and the 




4 5 6 7 8 9 

T(K) 

FigQ Experimental results for temperature depen- 
dence of the real part of AC susceptibility x'(T,h a c) 
for different AC field amplitudes h ac = 41.0, 59.6, 67.0, 
78.2 and 96.7mOe. 



measurements correspond to the zero-field cooling AC 
magnetic susceptibility. If we apply an AC external field 
H ac (t) = h ac cosujt normally to the 2D-JJA, then the 
total magnetic flux $(t) threading the four-junction su- 
perconducting loop is given by $(t) = $ ex t(t) + LI(t) 
where L is the loop inductance, & e xt{t) — SH ac {t) is the 
flux related to the applied magnetic field (with S ~ a 2 
being the projected area of the loop), and the circulat- 
ing current in the loop reads I(t) = Ic(T) sin </>(£). Here 
<p(t) is the gauge- invariant superconducting phase differ- 
ence across the ith junction. As is well-known, in the 
case of four junctions, the flux quantization condition 
reads Q31H3! 
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Fig|21 Theoretically predicted dependence of the nor- 
malized susceptibility x'(T,h ac )/xon (for better visual 
effect the curves are normalized differently with xoo = 
2-5xo, X03 = l-5xo, and xos = Xo where xo = S 2 /VL) 
on reduced temperature T/Tc according to Eqs.(3)-(5) 
for / = 0.5 and for "quantized" values of /3l(0) = 
2n(n + 1) (from top to bottom): n = 0, 3 and 5. 



where n is an integer and for simplicity we assume as 
usual that [2m] 4>i = 4>2 = 03 = 4>± = <P- 

To properly treat the magnetic properties of the sys- 
tem, let us introduce the following Hamiltonian 



H(t) = J(T)[l 



cos <f>(t)} + ^LI 2 (t) (2) 



which describes the tunneling (first term) and induc- 
tive (second term) contributions to the total energy of 
a single plaquette. Here, J(T) = (§ Q /2n)I c {T) is the 
Josephson coupling energy. 

Before turning to the interpretation of the observed 
step-like structure of x' (T, h ac ) , we would like to briefly 
comment on the origin of reentrant behavior in our un- 
shunted arrays (which has been discussed in much de- 
tail earlier, see El)- A comparative study of the 
magnetic properties of two-dimensional arrays of both 
unshunted and shunted (using a molybdenum shunt re- 
sistor short-circuiting each junction) Nb — AlO x — Nb 
junctions revealed |12j that the dynamic reentrance 
phenomenon can be explained by properly addressing 
the (neglected in our overdamped model) damping ef- 
fects associated with the finite value of the Stewart- 
McCumber parameter /3 C (T) = 2TtC J R 2 ,I c {T)/<S>n 
(where Cj is the capacitance and Rj is the quasi- 
particle resistance of the array). More precisely, the 
reentrance was found to take place in considered here 
unshunted arrays (with /3c(4.2if) = 30) and totally dis- 
appeared in shunted arrays (with /3c(4.2A") = 1). It is 



important to mention that both arrays had the same 
value of the (3l parameter, namely ^l(4.2K) = 30. 

Returning to the discussion of the observed here geo- 
metrical effects, we notice that the number of observed 
steps n (in our case n = 4) clearly hints at a possi- 
ble connection between the observed here phenomenon 
and flux quantization condition within a single four- 
junction plaquette. Indeed, the circulating in the loop 
current I(t) = Ic(T) sin <f>(t) passes through its maxi- 
mum value whenever <j>(t) reaches the value of ?(2n+ 1) 
with n = 0, 1, 2... As a result, the maximum number of 
fluxons threading a single plaquette (see Eq.(l)) over the 
period 2-k/u becomes equal to < >= (n + l)$o. In 
turn, the latter equation is equivalent to the following 
condition /3l(T) = 27r(n+ 1). Since this formula is valid 
for any temperature, we can rewrite it as a geometrical 
"quantization" condition /3_l(0) = 27r(n+l). Recall that 
in our array /3l(0) = 31.6 (extrapolated from its exper- 
imental value (3l(A.2K) = 30) which is a perfect match 
for the above " quantization" condition predicting n = 4 
for the number of steps in a single plaquette, in excellent 
agreement with the observations. 

Based on the above discussion, we conclude that in 
order to reproduce the observed temperature steps in 
the behavior of AC susceptibility, we need a particular 
solution to Eq.(l) for the phase difference in the form 
of 4> n (t) = |(2n + 1) + 5<t>(t) assuming 5(j)(t) < 1. Af- 
ter substituting this Ansatz into Eq.(l), we find that 
4> n {t) ~ fn+i/3 L (T) + i/cos^ where / = 2irSh ac /% 
is the AC field related frustration parameter. Using this 
effective phase difference, we can calculate the AC re- 
sponse of a single plaquette. Namely, the real part of 
susceptibility reads 



X'(T, h ac ) = - f d(ut) CO S (cJ%n(*) 
7T JO 



where 
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Here V is the sample's volume. 

For the explicit temperature dependence of fit, (T) = 
2irLIc(T)/§o we use the well-known |14II15| analytical 
approximation of the BCS gap parameter (valid for all 



temperatures), A(T) = A(0)tanh ^2.2^^ J with 

A(0) = l.76kBTc which governs the temperature de- 
pendence of the Josephson critical current 



Ic(T) = I C (0) 



A(T) 



A(0) 



tanh 
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(5) 



Fig. |21 depicts the predicted by Eqs.(3)-(5) dependence 
of the AC susceptibility on reduced temperature for 
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FigEl Fits (solid lines) of the experimental data for 
h ac = 41.0, 59.6, 67.0, 78.2, and 96.7mOe according to 
Eqs.(3)-(5) with I3 L (0) = 10tt. 

f = 0.5 and for different "quantized" values of /?l(0) = 
2ir(n+ 1). Notice the appearance of three and five steps 
for n — 3 and n = 5, respectively (as expected, case 
n = corresponds to a smooth temperature behavior 
without steps). 

In Fig. |3| we present fits (shown by solid lines) 
of the observed temperature dependence of the nor- 
malized susceptibility x'{T,h ac )/xo f° r different mag- 
netic fields h ac according to Eqs.(3)-(5) using /?l(0) = 
107T. As is seen, our simplified model based on a 
single-plaquette approximation demonstrates an excel- 
lent agreement with the observations. 

In summary, a step-like structure (accompanied 
by previously seen low-temperature reentrance phe- 
nomenon) has been observed for the first time in the 
temperature dependence of AC susceptibility in artifi- 
cially prepared two-dimensional Josephson Junction Ar- 
rays of unshunted Nb — AlO x — Nb junctions. The steps 
are shown to occur in arrays with the inductance related 
parameter Pl(T) matching the "quantization" condi- 
tion /3l(0) = 2n(n + 1) where n is the number of steps. 
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